Scaling and isospin effects in quasielastic lepton-nucleus scattering in the Relativistic 

Mean Field Approach 
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The role of isospin in quasielastic electron scattering and charge-changing neutrino reactions is 
investigated in the relativistic impulse approximation. We analyze proton and neutron scaling 
functions making use of various theoretical descriptions for the final-state interactions, focusing on 
the effects introduced by the presence of strong scalar and vector terms in the relativistic mean field 
approach. An explanation for the differences observed in the scaling functions evaluated from (e, e') 
and (v, fi) reactions is provided by invoking the differences in isoscalar and isovector contributions. 
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Extensive analyses of quasielastic (QE) inclusive elec- 
tron scattering data performed in recent years P, 0, H, 0| 
have clearly demonstrated the quality of the behavior 
known as scaling. These analyses are based on the so- 
called superscaling function, f(ijj'), obtained by divid- 
ing the cross section by an appropriate function which 
contains the single nucleon physics, and plotting the re- 
sult against the scaling variable ip'(q,u>) (see, e.g., [5J). 
One then studies the dependences upon the momentum 
transfer q and the specific nucleus chosen. From (e, e') 
world data one concludes that scaling of the first kind 
(no dependence on q) is reasonably respected at energies 
below the QE peak, whereas scaling of the second kind 
(no dependence on the nuclear species) is fulfilled very 
well in the same region. The simultaneous occurrence of 
both kinds of scaling is called superscaling. At energies 
above the QE peak, breaking of both kinds of scaling is 
observed, residing mostly in the transverse channel, and 
likely due to effects beyond the impulse approximation. 

Experimental data lead to a scaling function with a 
characteristic asymmetric shape, having a tail that ex- 
tends to high values of the transferred energy to (posi- 
tive values of the scaling variable tp'(q,ui)). The asym- 
metric shape of the scaling function is largely absent in 
non-relativistic (NR) models based on a mean field ap- 
proach. In contrast, the study presented in 0, 0] has 
shown that the asymmetry can in fact be obtained within 
the relativistic impulse approximation (RIA), given that 
a description of final-state interactions (FSI) using strong 
relativistic mean field (RMF) potentials is assumed. Re- 
cently, we have shown [8| that an asymmetrical scaling 
function can be also obtained within the framework of 
a semi-relativistic (SR) model, based on improved NR 
expansions of the on-shell electromagnetic current, pro- 
vided that the FSI are described by the Dirac equation- 



based (DEB) potential [9] derived from the RMF. Note 
that, in the SR model, the nonlocalities arising from the 
NR reduction of the Dirac equation — and incorporated 
into the wave functions through the Darwin factor - 
are essential to reproduce the asymmetric shape of the 
scaling function. 

The data analysis of the separated longitudinal (L) and 
transverse (T) contributions to the scaling function car- 
ried out in 0, 0| has shown that, whereas the L response 
does scale to a universal curve, the T strength increases 
with the transfer momentum q and/or the mass number 
A. This excess of strength in the transverse channel is 
not entirely understood, although different effects rang- 
ing from FSI effects to MEC contributions have been in- 
voked to explain it. This result also connects with the 
breaking of the zeroth-kind scaling (defined as the equal- 
ity of the scaling functions obtained from the separated 
L and T contributions, viz., = /t(V ; '1_ = fW)) 

observed in the data. From previous studies @, i, 03 it 
turned out that the zeroth-kind scaling is closely fulfilled 
by various models based on the impulse approximation. 
This is the case of the Relativistic Fermi Gas (RFG), 
by construction. The use of traditional NR and SR ap- 
proaches also leads to similar longitudinal and transverse 
scaling functions. This occurs for different descriptions 
of the FSI, namely, using the same Woods-Saxon poten- 
tial as in the initial state, which leads to symmetrical 
scaling functions that do not agree with experiment [l(| , 
and making use of the DEB potential that produces the 
correct amount of asymmetry in the scaling function [8[. 
Finally, the RIA also leads to the fulfilment of zeroth- 
kind scaling when the plane wave limit is assumed [7|, 
i.e., the final nucleon state is described as a free rela- 
tivistic on-shell particle. This is known as Relativistic 
Plane Wave Impulse Approximation (RPWIA). 
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Contrary to the above-mentioned models, breakdown 
of zeroth-kind scaling is observed in the RIA model with 
FSI described by the strong scalar and vector potentials 
of the RMF, which yields a scaling function found to be 
in excellent agreement with the data. The amount of 
scaling violation depends on the particular prescription 
chosen for the current operator. This is illustrated in the 
upper panels of Fig. [TJ where we show /(?/>') as well as 
its longitudinal fL(ip') and transverse fxiip') contribu- 
tions for 12 C. The momentum transfer has been fixed to 
q = 1 GeV/c when computing Jl.t, while, for the global 
scaling function /, in addition the beam energy has been 
fixed to e — 1 GeV, implying that the scattering angle is 
a function of tp' . The left and right panels correspond to 
results for the CC2 and CC1 prescriptions of the electro- 
magnetic current operator, respectively 0, [Hj]. I* 1 both 
cases breaking of zeroth-kind scaling is clearly observed, 
fx being about 20% larger than /j, for CC2 and almost 
twice as large for CC1. This result, first observed in 0, 
differs from the RPWIA and SR models considered in 
the literature 0, H, [H| where zeroth-kind scaling is well 
obeyed. Moreover, it turns out in the so-called scaling 
region (i// < 0) that, within the RMF model, scaling of 
the first kind is fulfilled separately by /l and fx with 
only modest scaling violations in each case. However, 
the two scaling functions obtained are different, with fx 
lying higher than //,. 
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FIG. 1: (Color online) Scaling function f(ip') compared with 
separate L and T contributions. All results correspond to the 
analysis of QE (e, e') on C. For the separate scaling func- 
tions the momentum transfer is q = 1 GeV/c, while for the 
(e, e') results the incident electron energy is s = 1 GeV. The 
results correspond to the RMF description of FSI (top pan- 
els) and the EM A approach (bottom). The current operators 
CC2 (left) and CC1 (right) have been considered. 



We should point out once more that, in contrast, the 
SR approach with a convenient description of FSI based 
on the DEB potential (derived from the same relativis- 
tic Hartree potential) yields results Q consistent with 



the asymmetry shown by the RMF model and data, but 
which respect zeroth-kind scaling almost exactly. Al- 
though extensive tests showing the reliability of the SR 
expansion [l2j have been performed within the context of 
the RFG and RPWIA approaches, the effects introduced 
by the dynamical enhancement of the lower components 
in the Dirac spinors (rJ Q [l5| , accounted for within a 
fully relativistic calculation, are not present in the SR 
approaches and this may be the reason that the fully 
relativistic and the SR calculations for the scaling func- 
tion differ. From previous studies 0,0]) we have shown 
that the scaling function is very slightly modified by the 
dynamical relativistic effects linked to the initial bound 
nucleon states. 

The presence of strong scalar and vector potentials 
in the nuclear states (mainly the final one) leads to a 
significant enhancement of the lower components of the 
four-spinors describing the nucleon wave functions in the 
relativistic approach. Due to this, the ratio between up- 
per and lower components of the fully relativistic calcu- 
lation is quite different from the one for free spinors or 
implied in the NR or SR approaches [HI, [l7| • In order 
to assess the influence of the enhancement of the lower 
components (also called spinor distortion 18]) the fully 
relativistic calculation is compared here with the effective 
momentum approach (EMA) [|, Q]| G3, [3 . In the EMA 
the relationship between upper and lower components 
is forced to be the same as for free spinors. Thus the 
EMA wave functions lack the dynamical enhancement of 
the lower components due to the presence of strong po- 
tentials. This is also the case of SR approaches where 
projections over positive-energy states and truncated ex- 
pressions of the current operator are considered. Lacking 
this kind of distortion, the EMA results should lie closer 
to the so-called factorized result 

EE 03. 

Scaling of the zeroth kind is well fullfilled in the EMA 
model, as is clearly illustrated in the bottom panels of 
Fig. [T] where we present the scaling functions again cor- 
responding to the two currents operators, now evalu- 
ated with the EMA. These results illustrate the crucial 
role played by spinor distortion in leading to zeroth-kind 
scaling violations. In fact, the EMA approach leads to 
f(ip') w /z,(-0') « /t("0') using either of the current op- 
erators. 

This result is in accordance with the study performed 
in the context of the SR approach and the DEB descrip- 
tion of FSI [8]. To conclude, notice that projecting out 
the negative-energy components yields an enhancement 
(decrease) of the L (T) contributions. The magnitude of 
these effects being moderate (strong) for CC2 (CC1) cur- 
rent operators. This is in consonance with the known ten- 
dency of the CC1 operator to enhance the effects of spinor 
distortion in the electromagnetic observables, while CC2 
shows a more moderate dependence on the amount of 
distortion 0, [H, US • 

One important application of superscaling has been 
suggested in [l9j within the context of making realistic 
predictions for charge-changing (CC) neutrino-nucleus 



3 



differential cross sections which, for example, are of in- 
terest in neutrino oscillation experiments. The validity 
of the superscaling hypothesis, i.e., the existence of a 
universal scaling function in electroweak processes, con- 
stitutes an essential result which is supported by various 
theoretical studies @, 0, S, E3] and gave rise to recent 
applications to neutrino studies [2(| [H| . The universal 
character of the scaling function is the basis of the Su- 
perScaling Analysis (SuSA) introduced in Within 
SuSA, the experimental superscaling function extracted 
from the analysis of (e, e') world data is used to recon- 
struct CC neutrino-nucleus cross sections. However, it 
is important to point out that the extraction of the ex- 
perimental scaling function refers only to the analysis of 
the longitudinal function, /l, whereas, in contrast, (z/, fi) 
reactions are totally dominated by the purely transverse 
Tyy +Taa and T VA channels. Thus, one may question 
the validity of using extracted from electron scatter- 
ing data to predict (y, fj,) cross sections, which are dom- 
inated by transverse responses. This issue is particu- 
larly relevant within theoretical frameworks which lead 
to fh{^') 7^ friip'), i-e., violation of zeroth-kind scaling. 
In what follows we present a study of this issue within 
the RMF approach. Our aim is to answer the above ques- 
tion and clarify the degree to which the scaling hypothesis 
does or does not work. 

In Fig. [2] we present the scaling functions obtained from 
the calculation of CC neutrino and antineutrino-nucleus 
reaction cross sections. The kinematics correspond to 
neutrino (antineutrino) energy fixed to s = 1 GeV and 
lcpton scattering angle 6^ — 45°. In each graph we show 
the results obtained with the RMF approach applied for 
neutrinos (solid line) and antineutrinos (dashed) com- 
pared with EMA (in this case only the curve for neutri- 
nos is presented as the one for antineutrinos is very sim- 
ilar). The separate analysis of the transverse channels, 
T and T', is also shown in middle and bottom panels, 
respectively. The longitudinal contribution to inclusive 
CC neutrino-nucleus scattering is negligible and there- 
fore is not shown. The usual relativistic single-nucleon 
expression for the charged-current operator [Til [l9| has 
been employed. From results in Fig. [21 several basic con- 
clusions emerge. First, the scaling functions for both 
processes (u,^r) and (y, ) almost coincide. Second, 
it is verified that f{4>') ~ /t("0') ~ /t'(VO to a high 
degree, i.e., the separate transverse responses contribute 
similarly to the global scaling function. Notice that this 
is valid, not only for the EMA approach, but also within 
the fully RMF model. Third, the dynamical enhance- 
ment of the lower components in the Dirac wave func- 
tions leads only to a slight modification of the scaling 
functions. These result are in accord with ones observed 
for (e, e') reactions in the CC2 case (see Fig. Q}. 

For the same kinematics, the basic difference between 
electron and neutrino scattering is the nature of the ex- 
changed vector boson, a virtual photon probing the elec- 
tromagnetic current in electron scattering, and a W 
probing the weak current in CC neutrino-nucleus pro- 
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FIG. 2: (Color online) Scaling functions evaluated from 
(u,fi~) and (T>,fi + ) reactions. The results correspond to 
SuCv) — 1 GeV and # M = 45°. Results from relativistic calcu- 
lations using the RMF are compared with EMA (see text for 
details). Top, middle and bottom panels refer to the global 
scaling function, the Tvv + Taa contribution and the Ty A 
one, respectively. 



cesses. As a consequence, whereas for inclusive (e, e') 
processes all nucleons (protons and neutrons) in the nu- 
cleus contribute, in the case of CC neutrino (antineu- 
trino) reactions only the neutrons (protons) in the target 
contribute to the inclusive cross section. 

In order to compare (e, e') and (y,[/r) reactions, it is 
convenient to separate the contributions of protons and 
neutrons to the (e, e') scaling function in the form 
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where Gp and Gn are the single-nucleon functions (in- 
corporating both the longitudinal L and transverse T 
contributions) defined in [J, 0| for protons and neutrons, 
respectively, and G = Gp + Gn- The proton fp and 



4 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 








^ o 

*♦-» 






3; 



. CC2 f\ 


RMF 




protons 




neutrons 






.5-1-0.5 0.5 


1 1.5 2 2.5 C 


■ cci 


RMF . 


prot. 




v \ neut. 







1.5-1-0.5 0.5 1 1.5 2 2.5 3 



■ CC2 * 


EMA 




Sl protons 











.5-1-0.5 0.5 1 1.5 2 2.5 3 



FIG. 3: (Color online) Proton and neutron contributions to 
the scaling function obtained from the analysis of 12 C(e,e'). 
Same kinematics as in Fig. 1. The results correspond to the 
RMF approach with CC2 (top panel) and CCI (middle). The 
EMA case is presented in the bottom panel. 



after the inclusion of dynamical relativistic distortion of 
the spinors, due mainly to the presence of strong relativis- 
tic potentials in the final state. Thus, the effect of spinor 
distortion on the scaling function is isospin-dependent, 
making the separate scaling functions for protons and 
neutrons appreciably different, fp ^ /at, that is, violat- 
ing scaling of the third kind. This is connected to the 
breaking of the zeroth-kind scaling observed in Fig. [1] in 
the full relativistic calculation. As will be shown later, 
the isoscalar/isovector terms in both the L and T chan- 
nels in (e, e') reactions also yield significant differences 
when compared with the purely isovector contributions 
involved in CC neutrino- nucleus scattering processes. 
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FIG. 4: (Color online) Longitudinal and transverse scal- 
ing functions for (e, e') compared with f(ip') evaluated from 
(v : fi~) and (77, fj, + ). All results have been evaluated with the 
RMF approach using the CC2 current operator. The kinemat- 
ics selected correspond to fixed values of the incident lepton 
energy, e = 1 GeV, and momentum transfer, q = 0.7 GeV/c. 
The averaged experimental function extracted from longitu- 
dinal electron scattering data is also shown [j. 



neutron /jv scaling functions are shown in Fig. [3] for the 
same kinematics as in Fig. [T] The total scaling func- 
tion, obtained as an average of fp and /at with weights 
Gp/G and Gn/G, respectively, is also shown in the fig- 
ure with solid lines. We show results for the RMF ap- 
proach and for both current operators CC2 (top panel) 
and CCI (middle). For completeness, the results cor- 
responding to EMA are also shown in the bottom panel 
(here we only consider CC2, since results for CCI are ba- 
sically the same). As one observes, in EMA the proton 
and neutron scaling functions are almost identical i.e., 
fp = /at; such a result will be called scaling of the third 
kind or isospin scaling. A similar result also occurs for 
the NR and SR calculations. In contrast, protons and 
neutrons yield different scaling functions in the RMF ap- 
proach, with magnitude dependent on the choice of cur- 
rent operator. The results of Fig. [3] indicate that the 
proton/neutron balance in f(ip') is significantly modified 



The universal character of the scaling function and 
its validity for electromagnetic and weak interactions is 
further analyzed in Fig. 2] Here we directly compare 
the functions and /tWO obtained from (e, e') 

cross sections with the ones corresponding to (z^ M ,/i~) 
and (77^,^+) reactions. All results have been evaluated 
within the RIA and making use of the RMF potential 
to describe FSI. The prescription of the current oper- 
ator is CC2 for both electromagnetic and weak interac- 
tions 0| ■ The kinematics correspond to fixed values of the 
incident lepton (electron, neutrino, antineutrino) ener- 
gies, 1 GeV, and transferred momentum q = 0.7 GeV/c. 
Similar results are obtained if, instead of fixing the mo- 
mentum transfer, we select a specific value of the lepton 
(electron or muon) scattering angle. The averaged QE 
phenomenological function obtained from the analysis of 
(e, e') data [4j is also included in Fig. 0] As observed, the 
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theoretical curve for /l(^>') follows the behavior of the 
data very closely (referred only to the analysis of the lon- 
gitudinal scaling function), and this proves the capability 
of the RIA combined with the RMF potential to describe 
(e, e') data in the longitudinal channel satisfactorily. On 
the contrary, the transverse contribution fx{ip') overesti- 
mates the data by ~ 20% even in the region close to the 
maximum, if>' w 0. This result, arising from zeroth-kind 
scaling violation in the RMF approach, is not in conflict 

with (e, e') data that indeed leaves room for effects of this 

n r~ii 

type (see the general discussion in [3, |J|). 

Concerning the scaling function obtained for neutrino 
(and antineutrino) scattering reactions, one observes that 
it is much more in accordance with /l^') and hence 
with the electron scattering longitudinal data, than with 
friip')- Within the context of our model this outcome 
reinforces the validity of the general assumption implied 
by SuSA i.e., the use of the phenomenological scaling 
function (extracted from the analysis of longitudinal QE 
electron scattering data) to predict CC neutrino-nucleus 
cross sections. However, it is also striking that f(4>') for 
and reactions, which are totally dominated by the 
purely transverse (T, T 1 ) channels, coincides with the Jl 
function of (e, e') instead of fx, in contrast to what one 
might expect. 

In order to understand these results, let us start by dis- 
cussing some basic differences between (e, e') and (u, /i) 
reactions. In the former, the longitudinal and transverse 
channels contribute importantly (at least for some kine- 
matics), and in both responses isoscalar and isovector 
form factors enter. However, as shown in Fig. [3J the bal- 
ance between isovector and isoscalar contributions may 
change in a significant way due to the strong dynamical 
enhancement of the lower components in the outgoing 
nucleon Dirac wave functions obtained with the RMF 
potential. In contrast, only purely isovector form factors 
enter in CC neutrino-nucleus scattering. Hence, it is im- 
portant to evaluate isovector and isoscalar contributions 
in (e, e') reactions and their effects on the scaling func- 
tions. In what follows we investigate how the functions 
fx and fx obtained from (e, e') RMF calculations change 
when the isoscalar form factors are removed. Notice that 
proceeding in this way, we force the (e, e') to be purely 
isovector, similar to what occurs for (u, jj,). 

The results of our analysis are presented in Fig. [5] for 
the same kinematics as in Fig. 01 Again, we compare 
the scaling functions fx and fx for electrons with those 
of neutrinos and antineutrinos. Experimental (e, e') data 
are also included for reference. Top and bottom panels 
refer to different assumptions concerning the electromag- 
netic form factors entering in (e, e') reactions. First, the 
curves fx and fx in the top panel have been obtained 
assuming G^j = —G P M , i.e., the proton and neutron mag- 
netic form factors are simply set equal in size and oppo- 
site in sign. In this way, we remove the isoscalar contri- 
bution in the magnetic form factor which, consequently, 
becomes purely isovector. The proton and neutron elec- 
tric form factors are not modified. Therefore, results for 
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FIG. 5: (Color online) Same as Fig. 4, but with modified 
isoscalar/isovector contributions via the nucleon form factors 
in the (e,e') case. Top panel: G\j = —G V M , i.e., magnetic 
form factor purely isovector. Bottom panel: Gjj = G P E = 0, 
i.e., no convective terms. 



(e, e') in the top panel of Fig. \E\ reflect the scaling func- 
tions where the isoscalar contribution only enters through 
the electric content of the nucleons. It is important to 
point out that the strength of the transverse nuclear re- 
sponse function R T increases significantly when isoscalar 
contributions in the magnetic form factor are removed. 
However, concerning the scaling functions, the results in 
Fig. when compared with Fig. QJ show that the dis- 
crepancy between fx and fx gets smaller because a vis- 
ible decrease occurs for fx- In other words, removing 
the isoscalar contribution in Gm leads to a weaker viola- 
tion of the zeroth-kind scaling property (within the RMF 
context). 

In the bottom panel of Fig. [5] we show the results cor- 
responding to no convective terms, i.e., the electric form 
factors for protons and neutrons (in the electromagnetic 
sector) are forced to be zero. This is a very drastic as- 
sumption which leads to nuclear electromagnetic longitu- 
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dinal responses R L being very close to zero, having only 
relativistic-order contributions involving Gm p or Gmu- 
The reason to consider this non-convective limit comes 
from the effects introduced by the isoscalar/isovector 
contributions in the electric form factors of the nucleon. 
Obviously, the neglect of convective terms yields neither 
isoscalar nor isovector contributions. While being aware 
of the important differences introduced in the cross sec- 
tions due to the assumption G P E — = 0, the analysis 
of scaling functions, constructed by taking the proper ra- 
tios between the nuclear responses and the single-nucleon 
ones, requires the use of the same approach (no isoscalar 
or no convective terms) in both the numerator (hadronic 
dependence) and denominator (single-nucleon). Hence, it 
is instructive to explore the behavior of the scaling func- 
tions in such approximations. The results in the bottom 
panel of Fig. [5] show that a unique (universal) scaling 
function emerges from the analysis of (e, e') calculated 
cross sections. Moreover, this function (without isoscalar 
terms) agrees with the one evaluated from {v^,^) and 
(puj/jr) processes and with f eX p(ip') extracted from the 
analysis of longitudinal (e, e') world data. 

In conclusion, we have investigated isospin effects in 
quasielastic electron and CC neutrino-nucleus cross sec- 
tions making use of the RIA and various descriptions of 
FSI with special emphasis placed on the RMF approach. 
Our main results can be summarized as follows: 

• The RMF description of FSI leads to a clear vio- 
lation of zeroth-kind scaling. This violation only 
occurs in this model, and not in other approaches 
based on non-relativistic, semi-relativistic or rela- 
tivistic plane-wave models. 

• The breaking of zeroth-kind scaling has been 
proven to be due to the important dynamical en- 
hancement of the lower components in the Dirac 
wave functions (mainly in the final state) produced 
by the use of strong scalar and vector relativistic 
potentials. 

• The balance between proton and neutron contri- 
butions in the scaling functions, namely third- 
kind scaling, evaluated from (e, e') reactions is also 
significantly affected by dynamical relativistic ef- 
fects. This contrasts with NR, SR and RPWIA ap- 
proaches, where protons and neutrons lead to sim- 
ilar scaling functions. The same comment applies 
to the effective momentum approach. 

• Finally, we have investigated in more depth the ori- 
gin of the differences observed between the scal- 



ing functions occurring in electron scattering and 
CC neutrino reactions. Contrary to what intu- 
ition would suggest, the longitudinal scaling func- 
tion for electron scattering is found to agree with 
the neutrino-scattering / (which is purely trans- 
verse) much better than does the transverse scaling 
function from electron scattering. We have shown 
that this result is consistent with the different roles 
played by isoscalar and isovector nucleon form fac- 
tors in the two processes. In part this result prob- 
ably arises because the convective effects are less 
important for the CC neutrino reaction than for 
electron scattering, especially for the axial-vector 
contributions in the former which are dominated 
by spin- flip matrix elements. 

These general results complement other previous find- 
ings @, 0, H3] and support the essential assumption of 
SuSA. Furthermore, the zeroth-kind scaling violations, 
not present in other models, may give us some clues as 
to how to proceed when trying to disentangle the sep- 
arate roles played by isoscalar and isovector form fac- 
tors. Obviously, more precision data, particularly for sep- 
arated longitudinal and transverse contributions to the 
cross sections and their effects on the scaling functions, 
are needed. The present experimental results indicate 
an excess of transverse strength below pion production 
threshold. It is important to identify the source of this 
strength (FSI effects in L and T channels, relativistic 
dynamics, effects beyond the impulse approximation...). 
The present study should be considered as a step in this 
direction. 
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